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We address the measurement of density of states within and beyond the superconducting gap in 
tunnel-coupled finite-size nanostructures using a capacitive method. Resonance shifts and harmonic 
generation are used to yield the full differential conductance spectrum without destruction of the 
low dimensionality otherwise induced by intimate ohmic coupling to an electrode. The method is 
particularly relevant to attempts to discern the presence of the fragile Majorana quasiparticle at 
the end of spin-orbit-coupled nanowires in appropriate magnetic field conditions by their signature 
mid-gap density of states. 
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INTRODUCTION 

The recent observation of a magnetic field-induced 
zero-bias conductance (ZBC) feature in the proximity su- 
perconductivity tunneling spectrum gap of a semiconduc- 
tor nanowire has spurred claims of a solid-state bound 
Majorana fermion. [l|-|3[ If true, this topological excita- 
tion may become the basis for realistic proposals of a 
fault-tolerant quantum computing scheme. [4|-[6f However, 
there still remain several intriguing inconsistencies in the 
available evidence before predictions [a, I7H13I] are unam- 
biguously confirmed. 

Perhaps most obvious is the failure of the supercon- 
ducting gap to close before the new ZBC feature ap- 
pears. This is particularly problematic because the gap 
closure marks a topologically critical boundary for the 
expected formation of the Majorana quasiparticle. Al- 
though recent theories have successfully argued that the 
discrepancy is caused by measurement of the wire-end 
local density of states (DOS) and not inconsistent with 
the topological phase, [14( a true measurement of "global" 
density of states gap collapse in magnetic field would pro- 
vide much-needed confirmation and significantly settle 
debate. 

Experimentally, this need to capture the global DOS of 
a one-dimensional nanowire presents a difficult measure- 
ment problem. Since the Majorana state is a property 
of the boundaries of the one-dimensional nanowire, we 
cannot simply contact it with a conductor to prepare a 
bulk tunnel conductance measurement and expect the 
dimensionality - and topological state - to persist. 

Here we propose a novel technique designed to solve 
this problem and measure the proximity superconduct- 
ing gap of a one-dimensional nanowire without any per- 
turbing ohmic contact. This approach, which exploits 
resonance and harmonic generation under AC excita- 
tion, is especially conducive to measuring the proper- 
ties of ensembles of quantum wires lithographically- or 
gate-defined from planar epitaxial heterostructures, and 
can be performed in a simple two-terminal geometry. 
To substantiate the proposal, a numerical simulation of 
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FIG. I. Side-view of the proposed hybrid device geometry, 
showing cross-section of nanowires capacitively-coupled from 
above, with tunnel conductance pathway to the substrate be- 
low. The equivalent circuit diagram, giving an equivalent 
capacitance C* , is overlayed on the left. 



harmonic generation using the discrete 1-d wire spec- 
trum as a function of magnetic field is performed. Fi- 
nally, a proof-of-principle experiment using a thin-film 
superconducting-normal metal tunnel junction is sug- 
gested to provide a means to benchmark against a trivial 
"hard gap" system. 



MEASUREMENT OF ZERO-BIAS 
CONDUCTANCE VIA RESONANCE 

The conductance of a tunnel junction (TJ) between 
a normal bulk 3-dimensional metal and a supercon- 
ductor reflects the superconducting gap in density of 
states near the Fermi energy, A.[15| In our case, the 
superconductor is a one-dimensional proximity-coupled 
nanowire and cannot be ohmically contacted with a 
DC probe. It can, however, be capacitively coupled to 
a counter electrode. This device, illustrated in side- 
view in Fig. 1, then has an equivalent circuit diagram 
where the tunnel contact to the nanowire consists of 
both a dissipative (G) and capacitive (Ctj) path in 
parallel. We imagine a lithographically-defined array 
of narrow-gap nanowires fabricated from a planar het- 
erostructure (such as Al/InAs/AlSb/GaSb as shown in 
the figure), patterned with Nb and selectively etched, 
encapsulated in conformally-deposited dielectric insula- 
tor (such as atomic-layer-deposited AI2O3), and capped 



2 




frequency [x 1 rad/sec] 



FIG. 2. Calculation of LC* circuit resonance shift when the 
zero-bias tunnel conductance to the wire varies for ui « G/C* . 
Panel (a) shows impedance magnitude, whereas panel (b) 
shows corresponding phase. Parameters chosen as L=100 /^H, 
d = 100 pF, Ctj = 400 pF, C s = 100 pF, and include a 
series resistance of 10 Q. Note that these parameters give 
\/L/C* ~ R (vertical dashed lines), so that cj ~ R/L. 



with a planar "top gate" metal counterelectrode capaci- 
tively coupled to the layers beneath it. We further note 
that giant bandgap bowing in the InAsSb system (l6j in- 
dicates that III-V materials with even larger spin-orbit 
interaction than InSb are available for investigation. 

When this system is driven in series with an inductor, 
its effective capacitance C* determines the resonant fre- 
quency wo = 1/V LC* , where L is the inductance. If L 
is chosen so that resonance occurs where the impedance 
of the TJ capacitive channel matches the zero-bias re- 
sistance of the dissipative channel R = 1/G, then the 
system is optimally sensitive to changes in zero-bias con- 
ductance G. For example, if the inductor is selected to 
match the normal state conditions at high temperature, 
then cooling the system below the critical temperature 
induces proximity effect from the superconductor, and 
the small-signal conductivity drops to a level determined 
by soft-gap creep, [ljj The current path then favors the 
capacitive channel in parallel, and the total capacitance 



C* drops because it is in series with the capacitor to the 
counter electrode. Experimentally, this is observed by a 
sudden shift of resonance to higher frequency. 

To illustrate this point, the magnitude of total 
impedance of an equivalent resonant circuit as a function 
of driving frequency for several resistance values is shown 
in Fig. 2(a). The corresponding phase is shown in Fig. 
2(b), where transitions from — ir/2 to ir/2 clearly iden- 
tify shifts in resonance to higher frequency as expected, 
approximately equal to the resonance width. 

This scheme can sense changes in zero-bias conduc- 
tance, for example, if and when a magnetic field Zee- 
man splitting induces Majorana formation at the center 
of the superconductor energy gap A. However, this zero- 
bias sensitivity is not enough: one needs to map out the 
closure of the gap even before the topological regime is 
reached. As we detail below, this can be determined by 
the harmonic generation obtained while driving the sys- 
tem in the nonlinear regime where voltage drop across 
the tunnel junction exceeds A/e. 



MEASUREMENT OF GAP CLOSURE VIA 
HARMONIC GENERATION 

Consider the simple case of a nonlinear resistor R and 
capacitor C in parallel, with an inductor in series. By 
summing the two parallel current paths and setting equal 
to series current through the inductor L, we have 



GiVi + G 3 V{ +iojCV 1 
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(1) 



where V± is the voltage across the parallel RC portion and 
V2 is the voltage across the inductor. By symmetry, only 
odd powers of voltage are present; G3 is the leading-order 
nonlinear component coefficient. By summing voltage, 
we also have 



V1+V2 = V(t) = V^l+iuLGi-ujZLC+iuiLGaV?). (2) 
If our driving voltage is 



V(t) = Ve iu \ 



(3) 



then our system will find steady state with a voltage 
across the RC device of 



Vi (t) = V^e*"* + V} 3) e i3ut (4) 
where the response at the driving frequency uj is 
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and the third harmonic amplitude is given by 



iojLG 3 (V^) 3 



jLGiV 3 



1 - lu 2 LC + itoLGi (1 - uj 2 LC + iuLGiY 

(6) 

Note that this response is linearly proportional to the 

nonlinear coefficient in the resistor conductance, G3. By 

(3) 

measuring the harmonic amplitude as a function of 
driving amplitude |V|, the higher-order nonlinear compo- 
nents present near the superconducting gap can therefore 
be determined. 



TIME-DOMAIN MODEL 

To model the relevant experimental aspects of this har- 
monic generation scheme, we have performed an explicit 
time-domain simulation of the high-frequency response 
of a simplified one-dimensional system incorporating the 
discrete spectrum of a finite wire. The geometry includes 
capacitive coupling from both the planar electrodes to 
the wire, and tunnel coupling to the grounded electrode, 
but ignores the direct electrode-electrode coupling. This 
enables the implementation of a one-dimensional finite- 
differences Poisson equation to algebraically determine 
the electrostatic potential V{t) on the floating wire: 



V R (t)-V(t) 

V 2 V{t) = - ° 



(d + D)/2 



e e(d + D)-A' 
(7) 

where D is the distance from the wire to the capacitive 
electrode at potential Vr{£) = Vsin(ujt), d is the distance 
to the tunnel electrode at ground, N(t) is the number of 
electrons charging the wire out of equilibrium, e is the 
dielectric permittivity between the wire and electrodes, 
and A is the cross-sectional area (wire length x width). 
This expression yields 
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(8) 



After this wire voltage is calculated, we can deter- 
mine the instantaneous current flowing into the wire 
I = V ■ G(V), where G(V) is the voltage-dependent con- 
ductance to the wire. In a timestep At, this current im- 
parts a transfer of charge and a change in N(t), requiring 
repeated calculation of wire voltage, then current flow, 
etc. at every timestep. 

From this wire voltage fluctuation in time, one can 
then determine the displacement current flowing between 
the capacitive electrode (at voltage Vr) and the wire (at 
voltage V) by using the Maxwell addition to Ampere's 
law 
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FIG. 3. Harmonic generation in time-domain simulation of 
the hybrid capacitor with nonlinear conductance spectrum at 
5Mrad/s as shown in the inset. Whereas the fundamental 
"W and second harmonic "2oj" signals are relatively unaf- 
fected, the 3oj response reflects the gap energy A, establishing 
the utility of this method to measure proximity gap closure 
of a nanowire in a magnetic field. 
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Comparison of the magnitude and phase of this signal 
to the sinusoidal driving voltage Vr{£) can be used to 
determine e.g. the complex impedance at a given simu- 
lation frequency u. The scheme has been benchmarked 
by considering a wire with continuous energy spectrum 
and linear current-voltage relationship with the tunnel 
contact, in which case the equivalent predictions of sim- 
ple circuit theory are asymptotically obtained in the limit 
At -> 0. 

We use the fast Fourier transform to determine the 
magnitude of third-harmonic 3w response of a one- 
dimensional wire. As discussed in the previous section, 
we expect this signal as a function of oscillation ampli- 
tude V at fundamental frequency uj to reflect the nonlin- 
earities in tunnel conductance spectrum. 

To test this approach, we simulate the response of the 
hybrid capacitor with a simple nonlinear gapped conduc- 
tance similar to what is expected from tunneling into a 
trivial superconductor. Here, the device geometry has 
d =10 nm, £>=100 nm, permittivity e = 12eo, and area= 
widthxlength = lOOnm x 1/im. The resistance varies 
from lGfl below the gap voltage to 0.1GJ7 above it. As 
shown in Fig. [3l no signature of underlying spectrum 
features can be seen in the fundamental "W response 
at 5MHz, and 2oj signals are absent as expected by sym- 
metry, even for multiple gap values A. However, the 
3w amplitude is nearly 1% of the fundamental, and ex- 
hibits a sharp threshold at a voltage amplitude which 
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FIG. 4. (a) Energy spectrum of the filled states of a spin- 
orbit-coupled nanowire relative to the Fermi energy as a func- 
tion of magnetic field, calculated by diagonalization of the 
finite-differences Bogoliubov-de Gennes hamiltonian. (b) Cal- 
culated amplitude of the third harmonic at 3oj in the time- 
dependent response from a capacitively coupled nanowire 
with this excitation spectrum, relative to fundamental at uj. 



scales linearly with A. We are therefore confident in 
interpreting the 3uj signal as a function of excitation am- 
plitude at uj to be proportional to the nonlinear com- 
ponents of tunnel conductance spectrum in this hybrid 
tunnel/capacitor geometry. 

The incorporation of a discrete nanowire spectrum into 
the time-domain simulation involves a modification of the 
tunnel current expression, which can be calculated easily 
at zero temperature: 



(10) 



where 7 is coupling energy from the tunnel electrode, 
is the Heaviside step function, and the sum counts all 
the Ei discrete energy states of the wire that are avail- 
able for transport due to occupation and Pauli exclu- 



sion. |\&f 1 2 is the electron- like probability of the wave- 
function with eigenenergy Ei, relevant for hamiltonians 
which include superconductivity pairing between elec- 
trons and holes. The excitation spectrum has been calcu- 
lated by numerical diagonalization of the single-band 1- 
dimensional Bogoliubov-deGennes hamiltonian [7, 18- 20j 
for a 1-micron-long wire with 40 spatial discretizations 
(yielding a 160x160 matrix encoded using sparse matrix 
manipulation in Matlab) with the following parameters: 
spin-orbit Rashba coefficient a =0.2 eVA, chemical po- 
tential jU =0 meV (at the spin-degeneracy Dirac point in 
the absence of magnetic field) , superconductivity pairing 
energy A =0.25 meV, effective mass m* = 0.026too, an d 
g-factor g = 50, appropriate for InSb. 

In Figure IU(a), we show the energy spectrum evolu- 
tion with magnetic field B. When the topological con- 
dition for a helical state g^BB > y fi 2 + A 2 is achieved, 
the superconducting gap disappears due to the coales- 
cence of energies from end-localized Majorana states; be- 
yond this magnetic field value, nearly degenerate states 
at the Fermi energy are seen. Due to the finite wire 
length, interactions between the states cause level repul- 
sion and oscillation - the so-called "smoking gun" of the 
Maiorana.plj 

Figure|U[b) shows the resulting 3w response of a hybrid 
capacitor using this spectrum and Eqn. [101 at 5 Mrad/s. 
Again, the device geometry has d =10 nm, D=100 nm, 
permittivity e = 12eo, and area= width x length = lOOnm 
x 1/mi. The tunnel coupling constant 7 is IneV, much 
smaller than the gap A, justifying our disregard for state 
broadening; this coupling value is equivalent to a tunnel 
resistance R «1 Gfl, which is necessarily large 

to constrain the flow of current in a single timestep of 
At =12.6 ns < 1 electron. 

All the spectrum features are captured directly in this 
3uj harmonic generation. The gap can be seen to close 
and a strong signal at low excitation amplitudes cor- 
responds to the creation of the Majorana at the Fermi 
energy. The strength of the third harmonic under con- 
ditions corresponding to the continuum/gap edge is ap- 
proximately 1% of the fundamental excitation amplitude. 
We expect this robust signal should be experimentally 
measurable with an appropriate lock-in amplifier or spec- 
trum analyzer. 



DISCUSSION 

The experimental feasibility of the approach to mea- 
sure tunnel conductance spectra outlined here can be es- 
tablished in a planar device with an equivalent electrical 
configuration, but with a trivial superconducting gap. 
For example, fabrication of a Cu/AlO^/Al tunnel junc- 
tion on an oxidized Si substrate with a thin insulating 
layer of Si02 will provide exactly the same configura- 
tion shown in Fig. 1, but upside-down: The capacitive 



counter electrode is now the substrate, and the tunnel 
contact is the top thin-film layer. 

This control experiment will be especially useful in re- 
fining the scheme, as independent contact to the pla- 
nar superconductor can be used to independently cali- 
brate using conventional conductance spectroscopy. At 
low temperatures, the zero-bias conductance as measured 
via the resonance method should drop considerably when 
cooling below the superconducting critical temperature, 
and the energy gap observed through harmonic genera- 
tion should close in modest magnetic fields, and no sig- 
nature of a Majorana ZBC should be seen. 

We also note that at sufficiently high DC bias, it may 
be possible to alter the chemical potential of the wire. 
This will be especially important to tune [i ~ to achieve 
the topological condition at magnetic fields low enough to 
preserve the superconducting state in the adjacent metal. 

The author acknowledges important conversations 
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generous hospitality of Prof. A. Yacoby at Harvard dur- 
ing the preparation of this manuscript. 
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